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ABSTRACT 


Molecular nanomachines are artificial or biological components which are in need to communicate with each other to perform a specific 
task. For this purpose, we consider two nanomachines called as transmitter nanomachine (TN) and receiver nanomachine (RN). TN 
transmits molecules in environment and RN receives the molecules through receptors, thus enabling them to communicate each other. 
The mathematical model which is based on ligand receptor binding is simulated in matlab code. The simulation results show that 
concentration of molecules, threshold concentration, pulse duration and temperature of environment affect on the channel capacity and 


thus appropriate values are be selected to achieve for maximum capacity of communication. B 
Keywords: Nanomachines, Molecular communication, Ligand receptor binding. 3] 
ne] 
1. INTRODUCTION 
Nanomachines or Nanosystems are objects with overall dimensions at or below the micrometer range and are made of assemblies of 
nanoscale components with individual dimensions ranging approximately between 1nm to 100nm. These nanosystems are in need to 
communicate/coordinate each other to perform a specific task. However, these nanodevices could not be communicated with the ov 
existing interconnection systems. Molecular communication is a recent interdisciplinary research area including nanotechnology, = 
biotechnology and communication technology using molecules as a communication carrier [Suda, et al., 2005, Albert et al. 1998]. In = & 
& 


nature, molecular communication is observed in living organisms to enable biological phenomena to communicate with each other 
[Lodish et al., 2000, Whitesides, 2001]. Hiyama et al explains the possibilities and research challenges for molecular communication 
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[Hiyama, et al., 2005]. In one paper, Suda et al introduce the concept of molecular communication and made an attempt for design of 
molecular communication system [Suda, et al., 2005, Nakano et al., 2007]. A communication system called gap junctions is described to 
communicate between two cells in [Akyildiz, et al., 2008]. In another technique, message bearing molecules in a small container known 
as vesicle conveys them along a filament connecting two devices using molecular motor [Hiyama, et al., 2005, Rospars, et al., 2000]. 
Cavalcauti et al propose a technique in which molecules may propagate in free space via Brownian motion sending the information in 
the pattern of molecule [Cavalcanti, et al., 2006]. There are reports on molecular communication based on ligand receptor bindings 
[Lansky, 2001, Moritani, Hiyama, Suda, 2006, Atakan, Akan, 2007, Lacasa, 2009]. In this paper, the main concern is to study molecular 
communication based on ligand — receptor interaction and find the optimum parameters such as concentration of emitted molecules, 
amplitude and duration of pulse of TN, total number of receptors on RN and temperature of the environment. The paper is organized as 
follows; in section 2, we briefly present the mathematical model available, the computational and simulation is discussed in section 3. 
The results of the simulation and interpretations are discussed in section 4. Conclusion of the overall work is given in section 5. 


2. MATHEMATICAL MODEL 


In ligand receptor binding model, the molecules L encounter receptors R and constitute complexes C (bound receptors) as, 
ky 
— 
fot Ie Z C 


If [L] and [R] are the ligand and receptor concentrations, then the rate of change of bound concentration [C] can be obviously written as, 
[Lansky,2001, Moritani, Hiyama, Suda, 2006, Atakan, Akan, 2007, Lacasa, 2009,] 


IC] _ 


= (k_, +k, [L])[C]+k,[L][N,] (1) 


where ki, k.1, are the binding rate and release rate respectively and [R]+[C] = [Na] is the total concentration of receptors available in 


RN. 


Now, equation (1) can be simplified as follows, 


k +k (Lyle k_ | +k IL]}t k_, +k [Lt 
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Considering the initial condition: t=O and [C(O)] =0 


k LIN, 
A = _ 

K+ KIL) 

k [LIIN, ] \k +k [ylt 
Ic) = ok it} 


kK +k,{Li| 


_lk 
Then, (C(p] = [C()]}1-e a (2) 


Here [ Co] is the steady state concentration of bound receptors, 


k LIN, | 
[Co] = —-—*— 
+k 


3. COMPUTATIONAL AND SIMULATION 
We assume TN emits molecules with concentration [L(t)] in the form of pulses with amplitude [L] (umol/liter) during to seconds with 


probability P, of releasing ligand molecule 
(L(t)] =[L] ; nt<t<ntot+t. (n=O, 1, 2, 3, ...... ) 
and [L(t)] = 0 with probability (1-P.) otherwise. (3) 


During the pulse time, concentration of [C(t)] (umol/liter) rises exponentially and at time tp when the pulse duration ends, [C(t)] starts to 


decay accordingly 


[C(t)] = [Coo] (1 - eK ,**,) forO<t< to 
[C(t)] = [Cto] e*.1"**) fort > to (4) 


In molecular communication paradigm, during the interval to, TN can emit either molecules (L) or transmit no molecule resulting in two 
molecular bits called MB1 and MBO. Consequently, if RN sense a concentration of molecules which is greater than a prescribed 
concentration [S] (umol/liter), RN decides TN has transmitted molecular bit MB1 during to. Conversely, if that concentration is less than 
[S], RN decides TN transmitted molecular bit MBO. The concentration of delivered molecules within to, [Ni] is given by 


IN, ]= f{C@dt 


tN, ]= [Ela 


el _ eID) de (5) 
ok, +k [L] 
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Since TN is continuously emitting every to, the previous delivered bits affect the concentration of molecules in the current interval. 
Hence, the concentration of complexes coming from the previous interval that still remain in the current interval can be given as follows 


[Np1=PIN,]fe*'dt © 
0 


Therefore, the expected total concentration of delivered molecules during to, will be [Nr]= [Ni] + [Nie]. Consequently the maximum 
probability of having success in transmission of a MB1 can be written as p1 = [Nz]/[S] and (1-p1) is the probability for receiving the 
erroneous MBO. Thus TN achieves to deliver MBO successfully with probability, p2 = [S]/[Nip] and it does deliver it incorrectly with 
probability (1-p2). It is also possible to detect erroneous molecular communication bits at RN side, by detecting MB1 when TN intended 
to transmit MBO or vice versa. Thus the molecule delivery capacity is defined as the maximum number of non-erroneous molecular bits 
which can be delivered within specific time duration. According to Pi, pi and p2, the channel can be modeled as a symmetric channel 


and then transition matrix M, is given by 


PP, (1-P,)(1-p,) 


Fel pi pAt=E) 
(7) 


Based on the transition matrix, the mutual information (MI) which states the number of distinguishable molecular bits, is as follows 


MI =-—| P, [Nr], q B-aS4 log| P, rl. a—p, {1 [S] } 
[S] [Nip] [S] [INP] 


-p,| al tog( (Nl) (1 Neel bog{ 1 7 Dal) 
[S] [S] [S] [S] 
[S] [S] [S] [S] 
1—P 1 1 1 1— 
Sos (anes) [ Ta oe( <a 


Finally, the capacity of molecular channel between TN and RN that is the maximum number of non erroneous molecular bits delivered 
within to, is max (MI) [Atakan, Akan, 2007, Lacasa, 2009]. 


(8) 


4. NUMERICAL RESULTS AND DISCUSSION 

For starting numerical calculation, we take the initial simulation parameter from the literature [Lansky,2001, Moritani, Hiyama, Suda, 
2006, Atakan, Akan, 2007, Lacasa, 2009]. First, we calculate MI varying with probability P. for different prescribe concentration [S] to fix 
appropriate concentration of molecules (threshold concentration) that must be delivered to RN within time interval to for a successful 
delivery of MB1. For this simulation, ki=0.1 pmol/liter/sec, [L]=1 umol/liter, [Nr]=0.001 pmol/liter, k-1+=0.0001 1/sec, to=1 sec are used. In 
figure 1, MI is shown with varying P, for different values of [S]=0.0001 umol/liter, 0.00001 umol/liter, 0.000001 pmol/liter. 
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Figure 1 
MI verses P, for different concentrations [S]. 


It shows that for lower concentrations [S]=0.00001 pmol/liter and 0.000001 pmol/liter, max(MI) increases with appropriate P\. It 
indicates that there is possibility to deliver non-erroneous MB1. However, for the higher concentration [S]=0.0001 ymol/liter, max(MI) is 
small and in the lower side of the probability Pi. From this analysis, we consider [S] in between 1.0*10°° umol/liter to 9.0*10°° umol/liter 
for further simulation. Figure 2 shows the plot of MI with variation of P, for different pulse amplitude [L]=0.5 umol/liter to 3 umol/liter. It 
is required to find appropriate amplitude such that TN can deliver sufficient concentration to RN. For this simulation we have taken [S] 
=1.0e-6 umol/liter. It is found that [L] smaller than 0.5 umol/liter cannot give maximum MI. However [L] greater than 0.5 umol/liter is 
sufficiently high to achieve maximum MI. Therefore, appropriate value of [L] must be selected to achieve maximum molecular 


communication capacity. 


a eee [LJ=0.5umol/liter 
--- [L]=1umol/liter 


r 3 
--:@&--- [L]=2umol/liter 
wv [LJ=3umol/liter 


Figure 2 
Ml verses P: for different pulse amplitude [L]. 


www.discoveryjournals.org | OPEN ACCESS 


pane | 9 3 


disc@very 


os the uns 


VOL.14, ISSUE 34, 2020 


DRUG DISCOVERY | ANALYSIS 


an. 
0.1 ws | ™ t=0.5sec --@-- t=1sec oe aw 
--M&-- t=1.5sec ---w- t=2sec gf} 
--@-- t=2.5sec --<4-- t\=3sec a 
0.0 . 1 1 1 1 1 1 1 1 
0.0 0.2 0.4 0.6 0.8 1.0 


P 


L 


Figure 3 
MI verse P for different pulse duration to. 
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Figure 4 
MI verses Pi for different concentration of receptors [Nr] on RN 


The pulse duration to is also an important simulation parameter. Figure 3 shows the plot of MI verses P. for different values of pulse 
duration to. For to=0.5 sec and 1 sec, maximum MI is obtained. However, to greater than 1sec, maximum MI decreases at the lower side 


of PL. It indicates that TN cannot deliver appropriate concentration of molecules to RN at higher pulse duration. 
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Figure 5 
MI verses P, for different temperature T of environment 


For concentrations of receptors on RN, [Nr]=0.01 pmol/liter to 0.00005 pmol/liter, the plot of MI with varying PLis shown in figure 4. 
It is seen that for lower values of [Ng]=0.00005 pmol/liter, and 0.0001 pmol/liter, the max(MI) is low at the higher probability. This is 
probable that although TN transmits sufficient number of molecules; there cannot be formed sufficient number of ligand receptor 
complexes on RN because of small concentration of receptors on RN. If we consider a higher concentration [Nr]=0.01 umol/liter, the 
maximum MI is still low. It means that TN may transfer more than the prescribe concentration to RN for delivery MBO results in an 
erroneous molecular MB1. In between the concentration in the range [Ne]=0.001 pmol/liter to 0.005 umol/liter shows maximum MI at 
appropriate probability. Therefore, it is required to select appropriate range of [Ne] to optimize communication capacity. Figure 5 shows 
the plot of MI with varying Pi for different temperatures T of environment. It is seen in the figure, for T= 300 K to 900 K, maximum MI 
can be achieved at an appropriate probability but as T increases, maximum MI decreases. It may be attributed due to increase of T, k-1 
decreases and ki increases such that TN can transmit higher concentration of molecules. On the other hand, TN cannot deliver 
concentration smaller than [S] for transmission of MB1. So here is the erroneous transfer of molecular bit. Hence, temperature of 
environment in which molecular transfer takes place is also an important parameter to achieve higher communication capacity. 


5. CONCLUSION 

In this paper, we have studied a simplified mathematical model for the possibility of molecular communication in the light of ligand 
receptor model. The simulation results revealed that concentration of emitted molecules, amplitude and duration of pulse of the 
transmitter nanomachine, total number of receptors on the receiver nanomachine and temperature of the environment affect the 
molecular communication capacity. We conclude that appropriate values of the above mentioned parameters are to be selected so that 


maximum channel capacity may be achieved. 
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